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ABSTRACT 

The thermal instability of a couple-stress fluid acted upon by 
uniform vertical rotation and heated from below is investigated. 
Following the linearized stability theory and normal mode 
analysis, the paper mathematically established the condition for 
Characterizing the nonoscillatory motions which may be neutral or 
unstable for rigid boundaries at the top and bottom of the fluid. It 
is established that all nondecaying slow motions starting from 
rest, in a couple-stress fluid of infinite horizontal extension and 
finite vertical depth, which is acted upon by uniform vertical 
rotation opposite to gravity and a constant vertical adverse 
temperature gradient, are necessarily nonoscillatory, in the regime 

where T A is the Taylor number and F is the couple-stress 

parameter, the result is also in accordance with corresponding 
configuration of Newtonian fluid when the couple-stress 
parameter F=0, Gupta et al. 4 
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1. INTRODUCTION 

Right from the conceptualizations of 
turbulence, instability of fluid flows is being 
regarded at its root. The thermal instability 
of a fluid layer with maintained adverse 
temperature gradient by heating the 
underside plays an important role in 



Geophysics, interiors of the Earth, 
Oceanography and Atmospheric Physics etc. 
A detailed account of the theoretical and 
experimental study of the onset of Benard 
Convection in Newtonian fluids, under 
varying assumptions of hydrodynamics and 
hydromagnetics, has been given by 
Chandrasekhar . The use of Boussinesq 
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approximation has been made throughout, 
which states that the density changes are 
disregarded in all other terms in the equation 
of motion except the external force term. 
Sharma et al. 7 has considered the effect of 
suspended particles on the onset of Benard 
convection in hydromagnetics. The fluid has 
been considered to be Newtonian in all 
above studies. With the growing importance 
of non-Newtonian fluids in modern 
technology and industries, the investigations 
on such fluids are desirable. Stokes 11 
proposed and postulated the theory of 
couple-stress fluid. One of the applications 
of couple-stress fluid is its use to the study 
of the mechanism of lubrication of synovial 
joints, which has become the object of 
scientific research. According to the theory 
of Stokes 11 , couple-stresses are found to 
appear in noticeable magnitude in fluids 
having very large molecules. Since the long 
chain hylauronic acid molecules are found as 
additives in synovial fluid, Walicki and 
Walicka 1 modeled synovial fluid as couple- 
stress fluid in human joints. An electrically 
conducting couple-stress fluid heated from 
below in porous medium in the presence of 
uniform horizontal magnetic field has been 
studied by Sharma and Sharma 10 . Sharma 
and Thakur 8 have studied the thermal 
convection in couple-stress fluid in porous 
medium in hydromagnetics. Sharma and 
Sharma 9 and Kumar and Kumar 3 have 
studied the effect of dust particles, magnetic 
field and rotation on couple-stress fluid 
heated from below and for the case of 
stationary convection, found that dust 
particles have destabilizing effect on the 
system, where as the rotation is found to 
have stabilizing effect on the system, 
however couple-stress and magnetic field 



are found to have both stabilizing and 
destabilizing effects under certain 
conditions. 

Keeping in mind the importance of 
non-Newtonian fluids, the present paper is 
an attempt to characterize the onset of 
instability analytically, in a layer of 
incompressible couple-stress fluid heated 
from below in the presence of uniform 
vertical rotation opposite to force field of 
gravity, when the bounding surfaces of 
infinite horizontal extension, at the top and 
bottom of the fluid are rigid. It is shown that 
for the configuration under consideration 



that, if 



4 + 2x 6 F + x*F 



< 1 , then an arbitrary 



neutral or unstable modes of the system are 
definitely nonoscillatory and, in particular 
the PES is valid, where T A is the Taylor 
number. 



2. FORMULATION OF THE PROBLEM 
AND PERTURBATION EQUATIONS 

Considered an infinite, horizontal, 
incompressible couple-stress fluid layer, of 
thickness d, heated from below so that, the 
temperature and density at the bottom 
surface z = 0 are T 0 , p Q respectively and at 
the upper surface z = d are T d , p d and that a 
uniform adverse temperature gradient 

M = — I is maintained. The fluid is acted 



dz 



upon by a uniform vertical rotation £2(0,0,£2). 

— > 

Let p , p, T and q(u,v,w) denote respectively 

the density, pressure, temperature and 
velocity of the fluid. Then the momentum 
balance, mass balance equations of the 
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couple-stress fluid (Stokes ; Chandrasekhar 
and Scanlon and Segel 5 ) are 



3 q 



q .V \q 



1 



Vp + g\l 



Pit 



Po 

V 2 q + 2] q x Q 



Sp_ 



(1) 



V.q=0 

The equation of state 
p = p Q [l-a{T-T 0 )] 



(2) 



Where the suffix zero refer to the values at 

— > 

the reference level z = 0. Here g(0,0,-g) is 
acceleration due to gravity. 

Let c v ,c pl denote the heat capacity 

of the fluid at constant volume and the heat 
capacity of the particles. Assuming that the 
particles and the fluid are in thermal 
equilibrium, the equation of heat conduction 
gives 



A) c v 

Or 

dT 

4 

dt 



dt 



q.V\T = q V l T 



q.V 

V J 



(4) 



The kinematic viscosity V, couple- 
stress viscosity// thermal diffusivity 

q 



K = - 



Po c v 



and coefficient of thermal 



Assume small perturbations around 
the basic solution and let Sp , Sp , d and 



q(u,v,w) denote respectively the 
perturbations in density, pressure p, 
temperature T and couple-stress fluid 
velocity (0,0,0). The change in density Sp 
caused mainly by the perturbation 9 in 
temperature is given by 
Sp = -ap o 0 . (6) 
Then the linearized perturbation 
(3) equations of the couple-stress fluid becomes 



d -f=-Usp-lae- 

dt p Q 









V 1 q+2\ qxQ. 


I A) J 








V.J = 0, 





(7) 



dt H 



(8) 
(9) 



Within the framework of Boussinesq 
approximation, equations (7) and (8), give 



3 V72 

— V w-ga 
dt 



v-^V 2 
Po 

dt dz 



d 2 e aV 

dx 2 dy 2 



V 4 w 



+ 2Q^ 

dz 



v-^Lv 2 

Po 



V 



(10) 

(11) 



expansion a are all assumed to be constants. 

The basic motionless solution is 
4=(0,0,0) , T = T 0 -jBz, H = (0,0,£l) 
and p = p 0 (l + afiz) . (5) denote the z-component of vorticity. 
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3. NORMAL MODE ANALYSIS 

Analyzing the disturbances into 
normal modes, we assume that the 
Perturbation quantities are of the form 



[w,0,g]=[w{z),®{z),z{z) 
Exp [ik x x + ik y y + nt) 



(12) 



Where k x ,k y are the wave numbers along 
the x and y-directions respectively 

k = (k x 2 +k y 2 p , is the resultant wave number 

and n is the growth rate which is, in general, 
a complex constant. 

Using (12), equations (9), (10) and (11), on 
using (8), in non-dimensional form, become 



(d 2 -a 2 ^tr + Fip 2 -a 2 f-ip 2 -a 2 ] 

w = _gad^e_^ dDZ 



l- F { D 2 -a 2 l D 2 -a 2 )-a\ = -JlA 



DW 



(d 2 - a 2 - Pl o)d ■■ 
where 



fid 2 



W 



(13) 

(14) 
(15) 



,, nd L v 

a = kd,a = , Pl = —,F- 

v K 



,T A = 



ia 2 d A 



,2 ' A 

Pod v 



Substituting w = W® 



z = 



© = ^-0 e and 

K 

^-Z e in equations (13), (14) and (15) 



and dropping (©) for convenience, in non- 
dimensional form becomes, 



w 



(D 2 -a 2 ^o + F{D 2 -a 2 f-(D 2 -a 2 )l 
= -Ra 2 &-T A DZ 

[{i-f(d 2 -a 2 )\D 2 -a 2 )-(r]z = -DW 
(D 2 -a 2 - Pl o)d = -W 



(16) 
(17) 

(18) 



Where R = ga ^ d , is the thermal Rayleigh 

KV 

number. 

Since both the boundaries rigid and 
are maintained at constant temperature, the 
perturbations in the temperature are zero at 
the boundaries. The appropriate boundary 
conditions with respect to which equations 
(16), (17) and (18) must be solved are 

W = DW = 0, 0 = 0 andZ = 0atz = 0 
andz = l. (19) 

Equations (16)-(18), along with 
boundary conditions (19), pose an 
eigenvalue problem for a and we wish to 
Characterize cr, when a r > 0 . 
We prove the following theorem: 



D = — and D @ = dD and dropping (©) for 

dz 



Theorem: If R) 0, F>0, T A }0, a r >0 and 
<r, * 0 then the necessary condition for the 
existence of non-trivial solution (w,&,z) of 
prandtl number, F is the couple-stress equations (16), (17) and (18) together with 



convenience. Here Pl =— , is the thermal 



parameter and T A is the Taylor number. 



boundary conditions (19) is that 
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\-f[d 2 -a 2 \D 2 -a 2 )-o*\z* =-DW* , (23) 

Therefore, using (23), we get 



Proof: Multiplying equation (16) by w* 

(the complex conjugate of W) throughout j^^^D^za^Jzl^-^-^-^f-^W, 



and integrating the resulting equation over o o o l > 

the vertical range of z, we get (24) 

r j , ,i \ j , ,ia 1 , Substituting (22) and (24) in the right hand 

cr \w*\D 2 -a 2 Wdz+F \w*\D 2 - a 2 ]Wd z - \W ., , V . & 

J 0 v ' J v ' J ^0) e ° ec l uatlon (20), we get 

[d 2 -a 2 Jwdz = -Ra 2 \w*®dz-T A \w*DZdz l , / , ,\ l / , , v? 

0 

Taking complex conjugate on both sides of w<fe - jw*^ 2 -a 2 f 
equation (18), we get 



o 

Wdz 



0 

,1 



(D 2 -a 2 - Pl <T*)®*=-W*, (21) 



Therefore, using (21), we get 



= Ra 2 \®[d 2 -a 2 - Pl (X*)~)*dz 

(25) 

-T A \z\{D 2 -a 2 )- F {D 2 -a 2 j-**\z*dz 



rTi7*<~.j (W™2 2 nil Integrating the terms on both sides of 

\W ®dz = -\®\D 2 -a 2 -p y CT p dz, (22) fe f 

J J equation (25) for an appropriate number of 

Also taking complex conjugate on both sides times b y makin § use of the appropriate 
of equation (17), we get boundary conditions (19), along with (17), 

we get 



\D 3 W\ +3a'\D z W\ +3a*\DW\ +a°\W\ \dz 



o\\\DW | 2 +a 2 |w| 2 }& +Fj- 
o L > 0 

+ j||dV| 2 + 2a 2 |DW | 2 + a 4 |w| 2 j<fe = Ra 2 }||Z3©| 2 + a 2 |0| 2 + Pl <7* |0| 2 jrfz - T A }{|DZ | 2 + a 2 |z| 2 }rfz 

-7- A Fj||D 2 z| 2 + 2 fl 2 |DZ| 2 + fl 4 |z| 2 U-r A <T*Jl Z | 2 * ^ 26 ) 



And equating imaginary parts on both sides We first note that since w and z satisfy 

of equation (26), and cancelling <r,(*0) w(0) = o=w(i) and z(0) = o = z(i) in addition to 

throughout from imaginary part, we get satisfying to governing equations and hence 

)\\Dw\ 2 + ay\ 2 \dz + Ra 2 Pl \\e\ 2 dz = T A )\z\ 2 dz , (27) we have from the Rayleigh-Ritz inequality 6 
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fpw| 2 &>;r 2 j|w| 2 &, 

0 0 

and 



i , l , 

l\DZ\ z dz>x 2 l\z\ z dz, 



1, , ,2 ,1 , 

(28) jp 2 w d z >7i 2 \\Dw\ d z and 



>l 9 I 2 7 1 , ,2 

jp 2 z &>^ 2 J|dz| z & , 



(30) 



(29) 



Further, 



for 



Further, multiplying equation (17) and its 
complex conjugate (23), and integrating by 
w(0) = o = w(i) and parts each term on both sides of the resulting 



Z(0) = 0 = Z(1), Banerjee et aC have equation for an appropriate number of times 



show that 



and making use of boundary condition on z 
namely z(0) = o = z(i) along with (17), we get 



1\\d 2 z\ 2 +2a 2 \DZ\ 2 +a 4 \zf \dz + F\ 



D 4 Z\ +4a 2 D 3 Z +6a 4 jD 2 Z +4a 6 ]\DZ\ +a S \z\ [<fe 
o' 1 0 



t 2F|||d 3 z| 2 +3a 2 |z) 2 z| 2 +3a 4 l\DZ\ 2 + a 6 \z\ 2 \dz + 2(J r ]\\DZ\ 2 + a 2 |z| 2 \d: 
.2\ n 2j 2 , 41^12 



0 L 0 J 0 

1 fi . ,2 „i „ i2 . ,1 „1 „ 1 



+ 2a r F\ \\D 4 z\ +2a 2 \D 2 z\ + a 4 |z| \dz + |a| J|z| <fe = \\dw\ dz 
0 I 1 1 1 J 0 0 



(31) 



Further, by utilizing boundary conditions 
(19) and equation (17), it follows that 



1, ,2 

jp 2 z dz = Real part of 
o 1 



-\DZ*D s Zdz 
0 



-\DZ*D i Zdz 
0 



\DZ*D i Zdz 
0 

< $\dz*d 3 zIi z , 
o' ' 

<}|z)ziz) 3 zU, 
0 I II I 

< jpz||D 3 zU , 

0 1 1 

1 



< 



\\DZ\dz 



l \ 7 I 2 
Jp 3 Z dz 

o 1 



_1 l 9 I 2 
j D 2 Z <fe 


2 


_1 l 1 I 2 
Jp 3 Z dz 


.0 




0 



(Utilizing inequality (30)), 
So that we have 



i 



'i^i 2 2 i h * i 2 

| D 2 Z dz\ <—\ UD 3 Z\ dz 
0 1 J T !n l 1 



Which yields 



h i 2 i'i i i 2 
jp 3 z Jz>^ 2 Jp 2 Z & , 



(32) 



Using inequality (29) and (30), inequality 
(32) becomes 



jp 3 z *>^ 6 ||z| z j z , 

O 1 1 0 

Also, 



1, ,2 

jp 3 z Real part of 



(33) 



-\D 2 Z*D A Zdz 
0 



(Utilizing Cauchy- Schwartz-inequality), 

Journal of Pure Applied and Industrial Physics Vol.1, Issue 4, 31 My, 2011, Pages (212-277) 



Ajaib S. Banyal, J. Pure Appl. & Ind. Phys. Vol.1 (4), 228-235 (2011) 



-\D 2 Z*D 4 Zdz 
0 

0 

<}\d 2 z*d 4 zIi z , 



< j\D 2 Z*\\D 4 Z\dz , 



do 2 zlD 4 Z 


dz, , 


o' " 


1 




"'l •> I 2 
l\D z Z\ dz 


2 


1, . ,2 
]\D 4 Z\ dz 


0 




o 1 



(Utilizing Cauchy- 
inequality), 

i 



Schwartz- 



2 I 



JjZ) j Z| cfej |j|D 4 Z| * 

(Utilizing inequality (32)), 
So that we have 



i_ 



\ \D 3 z\ dz 

Which yields 



\\D 4 Z\ dz 
O 1 



j \D 4 z\ dz>K 2 \\D i Z\ dz , 



(34) 



Using inequality (33), inequality (34) 
becomes 



\\D 4 z\ dz>7t r \\zfdz, 



(35) 



l\z\ 2 dz{ 



1 



\\Dwfdz , 



234 
(37) 



Now R > 0 and r A > 0, utilizing the 
inequalities (37), the equation (27) gives, 



l— i \ 

[x 4 + 2x 6 F + z & F}_ 



l\DW\ 2 dz + a 2 l\w\ 2 dz + Ra 2 Pl l\&\ 2 dz. 
0 0 0 



and therefore , we must have 



<o, 



(38) 



[ti 4 +2x 6 F + x % f) 

Hence, if 



;>!• 



a r >0 



and 



(39) 



then 



1>1 ■ 



-27t"F + x a F) 

And this completes the proof of the theorem. 

4. CONCLUSION 

This theorem mathematically 
established that the onset of instability in a 
couple-stress fluid in the presence of 
uniform vertical rotation, cannot manifest 
itself as oscillatory motions of growing 
amplitude if the Taylor number t a and the 
couple-stress parameter F, satisfy the 
inequality 



t <i. 



Now F ) 0 and ay > o , therefore the equation 
(31) gives, 



{ti 4 +2x 6 F + ti & f\ 

The essential content of the theorem, 
from the point of view of linear stability 
theory is that for the configuration of 
couple-stress fluid of infinite horizontal 
extension heated form below, having top and 
bottom bounding surfaces rigid, in the 
And on utilizing the inequalities (29), (30), presence of uniform vertical rotation parallel 



D 2 Z\ dz + 2F]\D 3 z\dz + F]\D 4 z\ dz(]\DW\ dz , (36) 
1 0 1 S 0 1 ' 0 



(33) and (35), inequality (36) gives 



to the force field of gravity, an arbitrary 
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neutral or unstable modes of the system are 
definitely nonoscillatory in character and in 

particular PES is valid, if-*— ? — v< l , 

the result is also in accordance with 
corresponding configuration of Newtonian 
fluid when the couple-stress parameter F=0, 
Gupta et al. 4 . 
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